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Abstract 



We study the long time existence theory for a non local flow associated 
to a free boundary problem for a trapped non liquid drop. The drop has 
free boundary components on two horizontal plates and its free energy is 
anisotropic and axially symmetric. For axially symmetric initial surfaces 
with sufficiently large volume, we show that the flow exists for all time. 

Numerical simulations of the curvature flow are presented. 



1 Introduction 

The evolution of interfaces of structured materials is of interest in a wide 
range of disiplines related to materials science [l^. Structured materials 
such as crystals, polycrystals and liquid crystals have a surface energy 
which is anisotropic; their energy density depends on the direction of 
the surface at each point. Over the years, various methods have been 
developed to track these interfaces, including the phase-fleld and level 
set methods. Here, we will consider a particular free boundary problem 
utilizing the anisotropic mean curvature flow. 

The mean (isotropic) curvature flow with constrained volume was con- 
sidered in [9]. In relation to the free boundary problem considered here, 
the papers of Athanassenas [3], |4] are particularly relevant . Also, in the 
recent paper F5^ , volume preserving mean curvature flow in a Riemannian 
setting is studied. 

Volume constrained anisotropic mean curvature flow for hypersurfaces 
was considered in [T], [TS]. In these papers, the emphasis is on the evolu- 
tion of closed convex hypersurfaces. 

Consider an anisotropic surface energy which assigns to a sufficiently 
smooth surface with unit normal v the value 
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The function 7 : -i^ R 
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is assumed to satisfy a convexity condition: 



the surface 
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which is known as the Wulff shape . In this paper, it will be assumed that 

(Wl) W is a smooth, uniformly convex surface of revolution with vertical 
rotation axis. 

(W2) W is symmetric with respect to reflection through the horizontal 
plane z = 0. 

Although we will not assume it in general, the following condition will 
also enter into discussion. 

(W3) The generating curve of W has non-decreasing curvature (with re- 
spect to the inward pointing normal) as a function of arc length on 
{z > 0} as one moves in an upward direction. 

The condition (W3) will be referred to as the curvature condition. 

Because of (Wl), the Gauss map of W defines a bijection of W onto 
S^. Therefore quantities defined on W can be expressed unambiguously 
on the sphere. In particular, the principal curvatures of W with respect 
to the inward pointing normal, i — 1,2 along the latitudinal and 
longitudinal directions, are given respectively by 

1/M2 = 7 - , = (1 - ^ih" + 'i-/fJ'2 ■ 

The uniform convexity means that fii > holds. Moreover, the position 
vector ^1 on W can be expressed in terms of the position vector u on 
as 

(3) ^1 = —iy + "/'(i'3)Es. 

Given an oriented embedded surface X : E -> R'' with unit normal 
field V, we define the Cahn- Hoffman field as the composition ^ o i/. 

The first variation of energy defines the anisotropic mean curvature A by 

(4) 5J^=- [ ASX-fd'E+'f ((,xSX-dX). 

Some local expressions for the anisotropic mean curvature are: 

A — -tr&cesAdi' = -(V ■ D7 - 2H'y), A := d^i = (D^7 -I- 7l)|s2 . 

Here D, denote the gradient and Hessian operators acting on functions 
on S^. These can always be identified with tensors on E by parallel 
translation in R since TpE = r^(p)S'^. Another useful expression is 

A = -V-C, 

where V- denotes the surface divergence. 

The problem we consider here is to understand the evolution to equi- 
librium of a drop of material trapped between two horizontal planes Hi, 
i = 0,1 located at heights 2 = and z — h respectively. The surface 
energy is assumed to be of the form |l|. It will also be assumed that: 

• the initial surface is axially symmetric, 

• the generating curve of the initial surface is a graph over the rotation 
axis. 
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• throughout the evolution, the surfaces meets the planes Hi orthog- 
onally. 

• the initial surface extends smoothly to an infinite periodic surface 
with period h. 

Of course, the third assumption is a consequence of the fourth, but we 
include it for clarity. 

It is easy to see that the axial symmetry is preserved by the evolution 
considered below. It is also consistent with the form of the known mini- 
mizer for sufficiently large volumes and neutral wetting along the interface 
of E with the planes, TT. For small volumes, the drop must disconnect 
during the minimization process. 

It will be shown below that if the initial volume is sufficiently large, 
then the generating curves of the evolving surfaces are also graphs for all 
times for which the flow is defined. 

Although we expect that a similar analysis can be carried out in arbi- 
trary dimensions, we have decided to concentrate on the case of surfaces 
in because of its obvious physical significance. 



2 Preliminary results 



We will restrict our attention to axially symmetric surfaces with vertical 
rotation axis. In addition, the generating curve will always be a graph of 
the form r = r{z) and we write the surface sis X — {r{z)e^^ , z) where we 
have identified the first two coordinates on with the complex plane C. 
The outward pointing normal to the surface is i/ = (1 -|-r^)~^^^(e**, — r^). 
Consider the flow: 

(5) ^^{k~K)u, 
with the boundary condition 

(6) i/ ■ £3 = , on as . 



From (141, we obtain 



(7) dtz = (A - A) 1/3 . 

It is evident from this formula that, in general, z must depend on t, so, 
following [7], we write z = z{C,,t) where we choose C, to lie in the interval 
[0, h\. We then express the immersion as 

X = X(C,t) = (r(z(C,t),t),z(C,i)). 

We get 



r^Zi -h ri = (A- A) ^1-1^1 
Combining this with ([7|, we get, using that u>, — —rz/\/l + r'j, 



.2 ' 
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i.e. 



(8) n = (A-A)\/l + ri . 

The evolution for r is identical with the evolution for r under the flow 

{dtX)^ = (A-Ao)i/, 

where _L denotes the normal component. 

Short time existence for (18 1 is standard, see for example Theorem 



8.1.1 of [12j . By the equation (12|, (13l given below, both A and A are 



expressible in terms of r and its derivatives. 

The admissible variations for the variational problem are those which 
keep the surface between the planes, SX • £3 = on 9E , and fix the three 
dimensional volume enclosed within the surface: 

SX-udT, = 0. 

s 

The infinitesimal generator (A — A)u of the evolution ([5|, ([6| clearly sat- 
isfies these conditions. 

By (|6|, along 9E, the normal u is contained in the horizontal equator 
of S^. By (W2), y'^ua = 0) = holds and so by ([3|, C and Xt are parallel 
along 9E. we then get from the first variation formula Q, 

dtJ^lJ^t] = - / A(A - A) dS = - / (A - A)^ dS , 

so the fiow decreases the anisotropic energy. 

Remark When the wetting is not neutral, it appears to be problematic 
to construct an analogous flow. In this case, the boundary condition for 
the minimizer ^ • £3 = Ci 7^ on dT, n Hi, where d are non zero constants 
related to the coupling constants for the wetting energy. This boundary 
condition is incompatible with the fiow given by (|5| maintaining the drop 
between the planes . 

For the boundary value problem we are considering, the morphology of 
the minimizer depends the initial volume and whether or not the condition 
(W3) holds. If (W3) holds, it was shown in |11| that for volumes greater 
than or equal to a critical value Vb, all stable equilibria must be cylinders. 
Below this value, any minimizer must either disconnect or loose contact 
with at least one of the supporting planes. 

When (W3) does not hold, numerical simulations for a particular class 
of functionals, [2], show that anisotropic unduloids may occur as stable 
equilibria for a certain range of volumes. However, for large volumes, 
only cylinders occur and for sufficiently small volumes, there is no stable 
connected surface spanning the two supporting planes. 



We will next show that of the initial volume is sufficiently large, the 
surfaces will not pinch off, i.e. the generating curves are bounded away 
from the rotation axis. 

The following lemma uses calibrations to show a minimizing property 
of graphs with A = 0. It is well known and its proof is given for complete- 
ness. 
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Lemma 2.1. Let T, be a surface with zero anisotropic mean curvature 
which can he represented as a graph over a planar domain Q.. Let S be a 
piecewise smooth oriented surface which is contained in the cylinder x R 
and which has the same boundary as E. Then, 

< T[S] 

holds. 

Remark. In saying that E and S share the same boundary, we mean 
that E — 5 is the oriented boundary of an oriented 3-chain. 
Proof. Let ^ denote the Cahn-Hoffman field of E. The condition that E 
has zero anisotropic mean curvature is expressed V • ^ = 0. 

Because E is a graph, we can extend ^ to a field ^ on f2 x R by making 
^ constant on all vertical lines through points in fl. This field will satisfy 
V ■ ^, where V denotes the divergence operator on R'^. Then, by the 
Stokes' Theorem, we get 

-F[E] = / j{us)d-E= / C-i^sdE 

i-usds^ [ i-usds 

E Js 

-f{us) dS = T[S] , 
where the inequality follows from 

Corollary 2.1. Let T be an axially symmetric anisotropic surface energy 
and let C C If^ be a circle. Let S be any piecewise smooth compact surface 
bounded by C and which is contained in the cylinder over the disc bounded 
by C . Also, let D be the flat disc bounded by C . Then 

holds. 

Proof. This follows immediately from the previous lemma using the 
fact that the disc has zero anisotropic mean curvature. 

Proposition 2.1. Let Eo be an initial axially symmetric surface enclosing 
a volume V and intersecting the .supporting planes orthogonally. Assume 
that 

(9) .F[Eo]<7(e3)^ 

holds and that the flow ^ exists for all t G [0, T). Then, then no pinching 
occurs. In particular 

(10) r > Co 
holds for t G [0, T), where 



7(63)— - -^Po] =: TTCo . 
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Proof. Recall that is non increasing. Let C denote the cylinder 

between Ho and Di enclosing the same volume V[Eo] as Eq. Let rc be 
the radius of this cylinder. If p :— r((", ti) < cq for some t\ <T and some 
C,, we can find an annular part of the surface Et^ bounded by two circles 
of radii p and radius rt^ > rc- From this piece of Ef ^ , we form a piecewise 
smooth disc type surface by filling the circle of radius p with a disc. 

By the previous corollary, we obtain 

J^[Eo] + TTCo > H^ti] + vrp^ > 7(e3)7rrtj > '){e-i)-Krl = 7(63) '^^'J"^ , 
which gives a contradiction. 



Lemma 2.2. Assume the conditions (Wl), (W2) hold for W and that 
E is an axially symmetric surface intersecting the planes Yii orthogonally 
for which holds. Then there exists ci = ci(co) such that 

ci > |A! 

holds. If, in addition (W3) holds, then we have 

> A > ci > -00 . 

Proof. The idea of the proof is essentially taken from [3]. 

We first note that the curvature condition (W3) can be expressed 

(11) ly-id^^fii > 0, . 

For an axially symmetric surface whose generating curve is a graph, 
the anisotropic mean curvature is given by 



(12) A=^ + ^ = ^ 

and its average value is 

/>r(l+r,2)V2 dz 



(13) A: 



r(l + r2)i/2 dz 



Recall that /Xi are the principal curvatures of the Wulff shape W with 
respect to the inward pointing normal so << /li < 00 holds. It then 
follows easily that 

r,'^^r(l + r2)V2d^ r'^=ldz 

~ Jo r(l + r2)i/2rf2 J^r{l + rfy/2dz~ 

From the boundary condition, we have = on 9Eo. Note that 
rzz/{l +r1) = (arctanr2)2. We get 

^1 o rr^ 



— dE = 2-K — ^ dz 
E Ml Jo Mi(l + r2) 

— 2iY I — (arctanr^)^ 
Jo Ml 

= —2ti I — (arctan Tz) — 27r / dz{ — )r(arctanrz) dz , 
Jo Ml Jo Ml 
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Note that, using Proposition ( |2.1[ ), we obtain 
— 27r f,f — (arctanr^) dz 

> — > C3(C0) 



First assume tliat (W3) holds. Note that since 1^3 = —rz/\/l + r|, we can 
write 

-2./a.(^Marctan..)d. = -2. / 9., (^Marctan(^=^)) 

= -27r / (^).(arctan(^=^))d^3 



< 



by (111, since ^3 and arctan(f3(l — u^) ^''^) have the same sign. Also, 
since 8^3^.1/ ii\ is uniformly bounded, we obtain that 

> 'y,^ > C4(co) , 

and the result follows by combining this with the previous inequalities. 

It (W3) is not assumed to hold, then we easily obtain 
|(fei//ii)| dE < cr dz dO for a constant c while the integrand in the 
denominator is counded below by r. q.e.d 



3 Evolution equations 

Again, the governing evolution equation is 

(14) dtX = {A-R)u. 

For any smooth variation 5X = ijju + T of a surface X, the correspond- 
ing pointwise variation of the anisotropic mean curvature is 

(15) 5A = J[tp] + VA • T . 
Here J is the self=adjoint elliptic operator given by 

(16) J[u] = V ■ AVu + {Adv, dv)u . 
Using that the evolution of X is given by ([sjl , we obtain 

(17) Kt-V ■ AVK = {dv,Adv){K-K) . 

We recall from [TTj that the normal v satisfies the equation 

(18) V • AVvj + {Adv, dv)uj = -VA ■ , j = 1, 2, 3 . 



This is a consequence of (151 and the translation invariance of the func- 
tional. 
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To compute the evolution of the normal, we use that for 5X = ipv + T, 
one has Sv = — V?/) + dvT. Since Xt = (A — A);/, we get 

(19) ut = -VA 

We define the parabolic operator 

Lemma 3.1. Define lo := ^1 + r| = 1/ \/l — z^l . Then 



,-,r 1 — 2Vaj ■ ylVo; cj , , , , , 
(20) -PH = 1- {Adv,dv)uj. 



Proof. By combining ( |18[ ) and ( |19[ ), we obtain the vector equation 

= -VA . 

For the surface of revolution, we have 



= (yi - !/| ,vz) , 

where we have identified the space of the first two coordinates with the 
complex plane. It follows from ( 18 1 that 

(21) V • AV{^l~vl)e'% = -{dv, Adv){^l~vle''' - VA , 

where the last term denotes the projection of VA onto . 

The metric on E is dS'^ = (1 + rl)dz'^ + r'^dO'^ := ai (g) ai + 02 ® ct2 
If f{z) is a differentiable function, then V/ • V6'dE = df A *d6 = fzdz A 
*Q2/r = /zdz A {-uj/r)dz = 0, and therefore Vf-VO^ 0. Also A6I dS = 
d*d9 = d* (Q2/r) = d{~oj/r)dz = 0. Finally IVei^dS" = de* A *d6l = 
df A {-uj/r)dz = (l/r^) dS so |V6lp = 1/r^ 

Using these formulas to expand out the left hand side of ( 21 1 and using 
that only depends on z, we have 



V -A^i^l-vle^") 

= e'^V • ylV(\/l-K?) + ie^'^^^^^^^^UlVef + A9) 



Combining this with (21 1, we get 



(22) V ■ ylV(,^/l - i/g) = - (Adt/, di.) A^/ l-uS- e"'" VA . 
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From this, it is easy to obtain 



= V-AV 



-j2= - (1 - uir'v ■ Avi^i^i) 



+ (Adu, du)uj + e~'^a; VA - V(l - uly^ ■ AV x 1 - u?, 
^2r^ V 

--^^ + (Adv, du)uj + e"'*tj^VA + - Vcj • AVw 
U2r^ CO 



(*), 



From ([l9|, we get (\/r^^^)t = — e '^VA which imphes 



= (1 - K')~'e~'''VA = a;'e~'VA , 



Combining this with (*) yields (201 

We state the following well known Maximum Principle, |15) . 

Theorem 3.1. // the operator u V ■ A\/[u] is elliptic and h : M x 
[t\,t2\ — R satisfies 

{dt-V ■ AV)h < , 

holds, then 

ma.xh(-,t) < max/i(-,ti) . 

M M 

In the case we are considering, ellipticity of u i— >■ V • j4V[m] follows 
from the convexity condition Wl. 

Proposition 3.1. Assume that the initial surface satisfies the estimate 
Then the generating curve is a graph for all time t for which the 
evolution ^ exists. 

Proof. As in [3], we use the boundary condition ([bJi and the assump- 
tions on the initial curve to extend the surfaces generated by the flow to 
periodic surfaces. 

For any constant c, we obtain from (201 
(23) 

(log(w) - ct)t - V ■ ylV(log a;) - ct) = c = {A du, dv) . 

r^lJ.2 <^ 

By Proposition ( |2.1[ ), there is a constant co = co(ti) such that r > co 
holds for t € [t,ti]. Hence, for a suitable constant c, we have (log(aj) — 
ct)t — V ■ j4V(log(a;) — ct) < 0. It follows from the Maximum Principal, 
that 

(24) max (log(uj) — ct) < max \oeuj . 

0<z<d,0<t<ti 0<z<d,t=0 



From the definition of uj we have from ( 24 1 that there exist constants 
C2(Eo), a2(Eo) such that if the flow exists for t G [0,T), then 

(25) to < cae"^' : 



9 



for < t < T. 

Recalling the definition of uj, the result follows, q.e.d. 



Proposition 3.2. Assume the conditions (Wl) and (W2) hold. Assume 
that the flow ^ is defined for < t < T. Then the curvatures of the 
surfaces X{t) remain bounded, i.e. there exists a constant C3{T). with 

(26) kl{t) + k1{t) =: \\diyf{t) < caiT) ,Vt,0 < t < T 

where ki(t) denote the principal curvatures at time t. 

First note that k2 = —{ruj)~^. By Proposition ( |2.1[ | we have r > co(T) 
and C4(T) > oJ > 1 hy the previous proposition and the definition of oj. It 
follows that 1^2! < C5(r). Since i — 1,2, are uniformly bounded below 
and above, {0 < a < fii < 1/a for some a £ R), it is enough to show the 
existence of a bound 

(27) A' < C6(r) 



for < t < T and then ( 26 1 will follow 



We recall the standard formula: 

(28) P[fg] = fP[g] + gP[f] - 2V/ ■ AVg , 

for sufficiently smooth functions / and g. The endomorphism field A 
75^7 + 7id is self-adjoint so that the last term is symmetric in /, g. 



From the last equation (17 1 and (20 1 we obtain 



2^2 

(29) P[uj^] = ^ - 2{Adiy, dv)u)'^ - 6VajAVaj , 



(30) P[A^] = 2{Adu, diy) A{A - A) - 2VA • AVA . 

For a suitable function h{x) and B G R both of which we determine 
later, we get 

P[h{uj^)(A^ + B)] 

2 2 

= (A^ + B)h'{uP)i - 2{Adu,du)u? - QVujAVuj) 

- (A^ + B)h"{u?)Vuj'^ ■ AVcj^ 

+ h{u'^){2{Adp,dv)A{A - A) - 2VA • AVA) - 2h'{u)'^)Vu? ■ AVA^ 
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We take h{x) := e"^ for a constant a to be determined later, and get 

P[e'"^\A^ + B)] 

= e"'"' ( (A^ + B)a{ - 2{Adiy, dv)uj'^ - GVcjAVcj) 
-(A^ + B)aVu? ■ AVu? + 2{Adv, diy) A{A - A) 



= e""' ( ^"^^ + h2a(A2 + B) + 2A2 - 2AA1 {Adu, du) 

+ [-6a^(A^ + B)Vuj'^ ■ AVlj^ - 2VA ■ AVA ~ 2aSIu? ■ AVA^' 
= (**). 

We have 

-2a(A^ + B) + 2A^ - 2AA < -2a(A2 + B) + 3A^ + A^ , 
using 2xy < + . Because of Lemma (1131), we have the existence of a 



constant Ci(r) with A'^ < (ci(T))'^. Hence by choosing a, B >> 0, we get 
that the term in (**) which includes the factor {Adv, di/) is non positive. 

Next notice that since the tensor A is positive definite and self adjoint, 
at each fixed point p £ E, we have an inner product defined by {u, v) = 
Au ■ V, u,v £ TpE. We then have using 2\{u, v)\ < (u, u) + (v, v), 

2a|Vcj^ • AVA^I = 4a|A||Vcc;^ ■ AVA| 

= 4a|A||(Va;^ VA) 

< 4a^ A^ ( Va;^ , Vo;^ ) + ( V A, VA) 

= 4a^AVaj^ ■ AVuj'^ + VA • AVA . 

This means that both terms in (**) between the square brackets are non 
positive for suitable choices of the constants a and B and so we can con- 
clude that 

P[e^-\A- + B)]<e--''^^^^^, 



holds. Now recall that by Proposition (2.1 1, we have r > co(Eo) > 0. It 
then follows that 

P[e-*"e'"*''(A' + = e"*" (^-Me""' (A' + B) + P[e""' (A^ + B)] 

i 2 



< e"*" ( -Me""' (A' + B) 

< 0, 



i2a(A^ + B] 



for a sufficiently large constant M. Recalling from Proposition (3.11 
that to < ■y/r+~c|(r) , we find that ( 27 1 follows from the Maximum Prin- 
ciple, q.e.d. 
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4 Higher order regularity 



In this section we obtain bounds on higher order derivatives of the surface. 
In similar problems involving mean curvature flow, these are obtained from 
bounds on higher order derivatives of the second fundamental form \dv\^ . 
In our case, a nice evolution equation for this quantity is unavailable. 
However, since we are only working with surfaces of revolution, we can 
use (301 instead to obtain the desired bounds. 

Proposition 4.1. For positive constants am = OmiT), Tm = TmiT), 

there holds 

(31) 

P[1V'"A1'] < -a^lV^+^Af +r„,(l + IV^Al^) ,r„ = r^dV^Aj^) ,j < m 
for Q<t<T. 

Lemma 4.1. Let f ~ f{s) be a sufficiently smooth function. Then 

(32) |V'"/(s)<iVr/(s)|' + i5i, 

where Ei only depends on \\/-' f\, j < m — 1. 

Proof. The proof is straightforward and is left to the reader. 
Proof of Proposition ^TTp At a fixed time t, we let ui := ds, 002 ~ rdd 
and we write, using the summation convention, 

dA — AjCjj , VdA = KjkOjjLOk , ... 

Then at any time t, we can write 

(33) IV-^Al^ = |A,,....,„J^ n 9'''"^ 

where, although the coframe is fixed, all quantities depend on t. Through- 
out the evolution, all the surfaces are axially symmetric, so g" = l/gu. 
For a normal variation X(t) = X + tijjv + the first variation of the 
metric is dtQij ~ —2il}hij where hij are the coefficients of the second fun- 
damental form. Since all the surfaces are axially symmetric, gij = SijQa 
and hij — Sijhu for all t. Using this, we obtain from (33 1 

aiV-Aj^ = 4H(A-A)|V"A|^ + 2A,„...,„(atA,„...,„) g'"^" 

fc — 1. .m 

= 4//(A-A)|V"A|^ + 2(V'"A,V™J[A]) 
~ 2{V™A,V'"(V- AVA)) , 

where ~ means that the quantities are equal up to terms of orders less 
than or equal to m. In the given frame, for a function /(s), 

(34) V-AV/= — a.(^). 

ruj fii 

In particular, since A only depends on s, we get 

(35) dtlV^Af ~ — (VrA)(Vr+'A) + $,„(vrA)(vr+'A) , 
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for suitable functions $m which are bounded by lower order derivatives. 
Since IV^Aj^ only depends on s, applying (341 gives 



— 94 — {V'"A,V«V"'A)) 
ruj fii 

— (V'"A, V.V.V^A) + — |vr+'A|2 

+ — ( — )s(V™A,V.V™A). 
rid fii 



From the last equation, ( 35 1 and Lemma 



4.1 



we can write 



vr^'AI" + ^^{V^A, V.V^A) 



for a suitable function which depends only on lower order derivatives. 
The second term can be bounded 



2 —2 



|-ivr+'A|2 



-a^lvr+'Al^ 



By applying lemma (4.11, we obtain 
P[|V™A|' 

where for suitable e, we can take, using ( |25[ ), 

(36) (J„, — (max(^i)y'l + ci(r) 

This shows that ([3T| holds, q.e.d. 



Theorem 4.1. There exist constants CmiT) such that 
(37) 1 V^Al' < C™(T) ,Vt < T. 



Proof of Th. ipTTp. Denote by S{m) the statement (37 1. The state- 
ment 5'(0) is just ( |27[ ) which was shown above. We assume S{j) holds for 
all j < m. 

Let ^rn ■■= IV^Al^. By ([28|, ([29| and ([3l|, we find, for positive 
constants a and B 



Pie"" (*™ + B)] 



< e 



+ (1/„ + B)(a[— - 2{Adv,dv)u? - 6VcjAVcj] - aVo;^ ■ AVo;^) 
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We can estimate the last term above using 
2a|V*™ • AVu?\ < e^vE-^+i + e' 
We choose e < am{T) and choose B >> so that —GBViuAViu + 



£~ a ■ AVuj < holds. This is possible since cj is bounded by (25 1. 
Recalling that {(jj/r) is also bounded in any finite time interval, so we 
arrive at 

Pie^^'i^fm + B)] < Ce°"'(*„ + B) , 

for t < T and a suitable constant C = C{T). Here we are using the 
induction hypothesis since the constant r™ = r,„(|V-'Ap), j < m, have 
been absorbed into the costant C{T). It then follows that 

P[e-'^*e''"'(*™ + B)] < 

holds for < t < T, so the result follows from the Maximum Principle, 
q.e.d. 



Proposition 4.2. // (31) holds then for suitable constants, the principal 
curvatures satisfy 

(38) IV^fcil^Cc™, \V'^k2\<c^ 

hold. 

Proof. For a surface of revolution, the Codazzi equations reduce to 

{k2)s = (fci - ki^Tsjr . 

It follows by an easy induction argument that an upper bound for |V^fc2| 
and hence |V'"fc2!, can be obtained from upper bounds on |V-'fci|, i = 1, 2, 
2 < m. For this one needs ( |10[ ) and the fact that derivatives of r of order 
j have upper bounds which depend on derivatives of fci of order < j — 2 
since Vss = kiZs. 

For a surface of revolution, the anisotropic principal curvatures are 
the functions := kj/fij, j = 1,2. These are the eigenvalues of the 
differential : TE -s> TW. Note that A = Ai + A2. 

We have (fc-,)s = {\j)sfJ,j+^jl-i'j {i')3{—ki)zs from which it follows easily 
that 

\'^Tkj\ < (max/ij)|V!rAj| + E2 , 
\VT\,\ < (max(lM))|Vrfc.i+i?2, 
where E2, E'2 depends on derivatives of the kj's of order less than or equal 
to m — 1. 

Finally, we have from the definition of A, 

|V'"Ai| < IV^AI + |V'"A2 

< |V™A| + (max(l/M,))|Vrfc2| + £2 

< |V™A1 +E3 + E'2, 

where E3 (and E'2) only depends on |V-'fci|, j < m — 1. Using induction, 



(37 1 and (261, the result follows, q.e.d. 
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Theorem 4.2. Assume that the initial surface satisfies Then the 

flow exists for all time. 

Proof. The result follows by a standard argument. Assume the 
flow exists on a finite time interval [OjTmax). Because of the uniform 



estimates given by (261 and (381, the flow can also be extended smoothly 
to t — Tniax. Then, by applying the local existence result, the flow can be 
extended to [0, Tmax + e) for some e > 0. q.e.d. 



Remark We show convergence of the anisotropic mean curvature 
to its mean. Note that 

and therefore ^ 

H^o] > [ [ {A-Af dEt dt . 
Jo Jst 

In particular 

(A - A)'' dEt ^ , 

as t — >■ 00 . 

5 Numerical Results 

Based on the different descriptions of the evolving surface, different meth- 
ods can be used to numerically solve a surface evolution equation. These 
include parametric, level set, and phase field methods (see [6]). Each 
method has its own advantage and disadvantage. We choose the para- 
metric method, which basically is a front tracking method, namely, the 
surface is evolved and tracked according to the surface evolution equation 

Instead of using fully implicit schemes for the temporal discretization 
of equation (|8|, a semi-implicit backward Euler method is used. The 
key of time forwarding in the semi-implicit scheme is to approximate 
the nonlinear terms in the equation by using the previously computed 
approximated solution, while the linear terms still need to be solved im- 
plicitly. Therefore we can avoid solving systems of nonlinear equations 
(for example, using the Newton's method) at each time step and thus 
the computational cost can be reduced. On the other hand, the implicit 
feature will increase the stability of the scheme so the restrictions on time 
step sizes can be loosened. 

As for the spatial discretization, a second-order finite difference 
method can be used. Finite element method using piecewise linear func- 
tions can also be used, and if the mass is lumped, it will be equivalent 
to the finite difference method. However, based on our numerical exper- 
iments with the above two methods, we choose to present the method 
of using cubic spline approximations. The spline approximation not only 
provides a higher order method, it also ensures the continuity of the second 
order derivatives across the spatial nodes, which we think it is important 
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in the approximation of curvature flows. On the other hand, since the 
curvature flow we study here is essentially one dimensional, as we will see 
below, when using the spline approximation, the resulting linear system 
that we need to solve is very sparse and the computational cost will be of 
the same order as the cost using flnite difference or finite element method. 



5.1 Approximation by Splines 

In this section, we discuss numerical methods for simulating the ini- 
tial boundary value problem ([sjl and For a representative class of 
examples, we concentrate on the Rapini-Papoular functionals given by 
7 = 1 + evi , where — 1 < e < 1 . 

We use the clamped cubic spline approximation to numerically solve 
the evolution equation. Let T > and let TV, K be positive integers. Let 
= zi < ■ ■ ■ < Zff < Z]v+i = 1 be a partition of the interval [0, 1] and 
let = to < ■ ■ ■ < tx <= T be an equally spaced partition of the time 
interval [0,r]. 

We use the following notation 

r"^ = r{tk, z„), {rz)'iT, = rz{tk,Zn), (r^^)'^ = r^^{tk, z„), A'^ = A{tk, Zn), 

where < k < K and 1 < n < A'^ + 1. We also set 

r = T/K, Q{r(t,z)) = l + {r4t,z)f . 

Using the semi-implicit backward Euler method and the above nota- 
tions, the discretized evolution equation is 

for < fe < and 1 < < TV -I- 1. 

For an approximation of the generating curve of the axially surface 
r{t,z), we seek a spline function S{t,z) that satisfies equation (39 1 and 
the following properties : 

(i) S„{t,z) = ^(i,^)!^^ ^ J, the restriction of S{t,z) on interval 

[zn,Zn+i], 1 < n < TV, is a polynomial of degree no more than 
3. 

(ii) The first and second order partial derivatives of S with respect to z 
exist at nodes zi, - ■ ■ , zn+i and are continuous at the internal nodes 

Z2,--- ,Zn- 

(iii) = Q = ^(i,0) and —{t,ZN+i) = /3 = -^{t,l), where 
oz az oz az 

Q = /3 = if the contact angles of the surface at the top and bottom 

planes are required to be right angles. 



At time t — tk+i, < k < K ~ I, the unknown for equation (39 1 is 
denoted by a column vector of length 2TV: 

(40) x''+' = . . . , r%\\;d',+\- ■ ■ , 4+^)'' 
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where 



r^+^ = S(tfc+i, Zn), = l<n<iV + l. 

oz 



To derive a linear system of x^'^^ from equation (391, for each 1 < 
n < N, we adopt the following notations: 



h„ 

Sn 

s 



= Snitk, z) = S{tk, 2)| 

— Zn-\-\ 



{tk,z), 0<k<K, 



Z G [Zn+l,Z„] 



The piecewise-defined spline function S{tk, z) consists of the functions 
Sn of the following on the interval [zn, Zn+i]: 



(41) 



3h„s' - fc ^hi- 3h,y + 2s^ fc 



hi 



hi 



+ 



S {S- /in) ^fc , S(S - h„) 



hi 



n + l 



for < k < K . Following the standard theory about cubic splines, we 
derive the following linear system 



(42) 



h„ 



+ 



h„ 



+ l -I- h 

' n + l /i.„a„_]_ 



+ 2(/i„_i + /i„)d*+' + hn-^d^Xl =0. 2 < n < iV . 
Invoking the boundary conditions, the above system can be written as 
(43) Mix^bi, 



where x is given as in equation (40 1, 6i is a column vector of length A*' — 1 
given by 

h[ = [-/i2a,0,--- ,0,-/ijv/3]^ , 
and Afi = [Mn, A/12] is an (A^ - 1) x (2iV) matrix with 



A/ii 



and 



Mv 



hi 



V/l: 



2 3/M 



h3 



h2 

/l2 h2 



3 

/ljV-1 



A _ 3 

/iJV ftjV-1 



hN 



(]V-l)x{iV+l) , 



■2(/li+/l2) /ii 

h-s 2{h2 + /13) 



ftjV-l 2(/l]V_l + /iiv) hM-2 

hN 2(/ijv_i + /liv). 



(]V-l)x(]V-l) . 
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On the other hand, from the equation (39 1, we can derive the rest of 
the equations needed for solving x. For convenience, we use 



1 + irz{tk,z„)) , 



and 



(Ail)n Qn ' 



( ^ + A* 



for < k < K and 1 < n < A'^ + 1. Then, the system (391 can be written 
as 



(44) 



fc + l / \k-\-l /-k . k I k 



Since from equation (41 1, we have 



fc+i 



fc+i 



2d. 



fc+i 



Ad. 



fc+i 



-65 



fe+i 



hn. 



fe+i 



1 < n < iV + 1. 

1 < n < iV; 
n = N + 1, 



using the boundary conditions, we can write (|44| as 
(45) M2X = b2, 

where 62 is a column vector of length A'^ + 1 given by 

/ii(m' + ?-f)-4<f 

&2 = : 

/i]v-i(?)^-i +rJv_i) 

hM{n% + r%)-2peN 
./iiv(r;^+i+r{v+i) + 4/3C^+i. 

and M2 = [M21, M22] is an (A^ - 1) x (2A'') matrix with 



M21 = 



and 



hi + 



6cr 



6e 



6?^ 



6^^ 



JV+l 



flN + 



(iV+l)x(JV+l) , 



M2' 



4ci 2ei 



2CAr-l 

-2C: 



JV+l 



(iV+l)x(JV-l) . 



Combining systems (431 and (451, we finally obtain the linear system 
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(46) Mx = b, 

where 

^ [Mil M12I r ^ \bi 

[M21 M22J ' [&; ■ 

As we can see, the system ( 46 1 is very sparse and can be solved by stan- 



dard sparse solvers. We also would like to mention that the computation 

— k 

of A„ is done by using Gaussian quadratures. 

5.2 Numerical Experiments 

Some of the results from our numerical experiments will be presented 
below to illustrate the evolution the curvature flow that we have studied. 
In all of our simulations, we let z G [0, 1]. 

Example 5.1. In this experiment, e = 0.2, T = 3, iV = 500, r = 10~*. 
The initial profile of r is chosen to be the cubic Hermite interpolant that 
satisfies 

r(0,0) = 0.7, r(0, 1) = 0.4, r^(0, 0) = r;,(0, 1) = 0. 

The results are shown in Figure 1: 

(a) r{t, z), the profiles of the generating curve of the surface at different 
times are shown in plot (a). The initial profile is eventually evolves 
to a cylinder and remains thereafter. 

(b) The snapshots of the surfaces in three dimensional space at different 
times are also shown in plot (b). 

(c) The history of the values of the energy functional F (in equation 
|l|) is shown in plot (c), and it can be seen that the energy is 
decreasing as the surface evolves under equation ([sjl until it remains 
nearly unchanged, which numerically implies that a minimum of the 
energy has been reached and the minimizer is corresponding to the 
the surface of a cylinder. 

(d) The history of the values of the volume enclosed by the surface is 
shown in plot (d). Although it is seen that the volume is not pre- 
served at the beginning of the evolution process, the error (relative 
to the initial volume) is within 0.007%. We think this is due to the 
error in the numerical approximations. First, the volumes are com- 
puted by using Gaussian quadrature. Secondly, the quantity A in 

is approximated by A„ as in our semi-implicit scheme, 
which numerically violates the law of volume preserving curvature 
flow (equation ([8|), unless in the later stage of the evolution, the 
surface nearly has constant anisotropic mean curvature. 

According to Theorem 5.1 of [llj, the threshold of stability for cylin- 
ders is 

(47) 

^ ' p2(0) r^ - 
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For 7=1 + and e = 0.2, the cylinder is stable provided 



Vl + 2e _ 
r > ^ 0.3766. 

This experiment numerically verifies the above stability analysis. In other 
similar experiments with larger volume fractions, we have also observed 
that no pinching has occurred. 



Profiles, r(t,z). s=0.2, wO=0,w1=0 
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Figure 1: Exaniple |5.1[ Snapshots of surface profiles, values of energy functional, 
and volumes, e = 0.2, r(0, 0) = 0.7, r(0, 1) = 0.4. 

Example 5.2. In this experiment, e = 0.2, T = 3, iV = 500, r = 10"^. 
The initial profile of r is chosen to be the cubic Hermite interpolant that 
satisfies 



r(0,0) = 0.4, r(0, 1) = 0.2, r,(0, 0) = r,(0, 1) = 0. 
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The results are shown in Figure [2] The initial volume is about 0.2980. 
Comparing with 0.9847, the volume in Example |5.1[ this one is much 
smaller. However, the initial surface still evolves to a cylinder. Also, the 
errors in volumes (relative to the initial volume) are within 7 x 10"* %. 



Profiles, r(t,z). s=0.2, wO=0,w1=0 
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Figure 2: Exaniple |5.2[ Snapshots of surface profiles, values of energy functional, 
and volumes, e = 0.2, r(0, 0) = 0.4, r(0, 1) 0.2. 



Example 5.3. In this experiment, e = 0.2, T = 3, iV = 500, r = 10"^. 
The initial profile of r is chosen to be the cubic Hermite interpolant that 
satisfies 

r(0,0) = 0.3, r(0, 1) = 0.2, r,(0, 0) = r,(0, 1) = 0. 



The only difference between this example and Example |5.2| is the values 
of r(0,0). The results are shown in Figure [s] However, this time it can 
be seen that singularity occurs. 
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Example 5.4. In this experiment, e = 0.2, T = 4, iV = 1000, r = lO"*^. 
The initial profile of r is chosen to be the cubic Hermite interpolant that 
satisfies 

r(0,0) = 0.9, r(0, 1) = 0.1, r,(0, 0) = r,(0, 1) = 0. 
The result is shown in Figure [7j 

Example 5.5. In this experiment, e = 0.4, T ^ A, N = 500, r = 10~^. 
The initial profile of r is chosen to be the cubic Hermite interpolant that 
satisfies 

r(0,0) = 0.3, r(0,l) = 0.2, r,(0, 0) = r,(0, 1) = 0. 

The results are shown in Figure 

Through these examples |5.2| t o [5. 5[ we see that for small volume frac- 
tions, namely, when inequality ([9|) is not satisfied, the stability of cylinders 
is lost and singularity may develop, and therefore the flow may not exit 
for all time. 

We present two more examples: In Example |5.6[ e < 0, in Example 
|5.7[ a different initial profile is used. The results are shown in Figure [6] 
and Figure [Tj respectively. 

Example 5.6. In this experiment, e = -0.2, T = 2, TV = 500, r = 10"*. 
The initial profile of r is chosen to be the cubic Hermite interpolant that 
satisfies 

r(0,0) = 0.8, r(0, 1) = 0.3, r,(0, 0) = r,(0, 1) = 0. 
Example 5.7. In this experiment, e = 0.2, the initial profile is given by 
r(0, 2) = 1 + - cos(87r2:). 
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